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The multifractal scaling exponents are calculated for the critical wave function of a two- 
dimensional Dirac fermion in the presence of a random magnetic field. It is shown that the problem 
of calculating the multifractal spectrum maps into the thermodynamics of a static particle in a 
random potential. The multifractal exponents are simply given in terms of thermodynamic func- 
tions, such as free energy and entropy, which are argued to be self-averaging in the thermodynamic 
limit. These thermodynamic functions are shown to coincide exactly with those of a Generalized 
Random Energy Model, in agreement with previous results obtained using Gaussian field theories 
in an ultrametric space. 
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I. INTRODUCTION 

In recent years it has become clear that the wave func- 
tions of noninteracting disordered systems at a contin- 
uous metal-insulator transition have multifractal scaling 
properties. As opposed to a simple fractal, these sta- 
tistical self-similar wave functions cannot be described 
by a single fractal dimension, but instead an infinite set 
of scaling exponents is needed. Indeed, such families 
of scaling exponents have been obtained for the critical 
wave function at a localization-delocalization transition 
within several different frameworks, including perturba- 
tive renormalization-group treatments on replicated |l]] 
and supersymmetric 0j nonlinear sigma models, as well 
as numerical simulations |3| . Although these results seem 
to provide sufficient evidence for multifractality at the 
metal-insulator transition, none of them allows one to 
understand the full spectrum of multifractal exponents. 
This is so because nonperturbative techniques are needed 
to probe the full spectrum. 

Given the present state of affairs of this problem it 
would be highly desirable to have an exactly solvable 
system that exhibits a multifractal wave function. In 
recent years it has been shown that the Dirac equation 
in random fields in two-spatial dimensions (2D) is ac- 
tually an example of such a system. Moreover, it has 
also become clear that the random Dirac Hamiltonian in 
2D describes the universality class of the metal-insulator 



transition at a half-filled Landau level of disordered non- 
relativistic noninteracting electrons in a very high mag- 
netic field |@J|. This surprising result has been recently 
established by direct derivation || of the random Dirac 
equation from the Chalker-Coddington model of perco- 
lating edge states JjJ. 

Recently, exact results for the full multifractal spec- 
trum have been obtained for Dirac fermions interacting 
with a random magnetic field (§) . These results were ob- 
tained by mapping the problem to a Gaussian field theory 
in an ultrametric space (a Cayley tree). The calculation 
of multifractal scaling properties was then reduced to 
the computation of thermodynamic functions of a spe- 
cial generalized random energy model (GREM), which 
are known exactly [sHl^]. The field theory in the local- 
ization problem is defined in Euclidean space. Neverthe- 
less, it was argued that as the exponents are a measure 
of a global property, the difference between the ultramet- 
ric and the Euclidean metric would not alter the results. 
Moreover, a phase transition that occurred in the ultra- 
metric model 1 11 also seems to be present in the original 
Euclidean problem, as evidenced by analytical arguments 
and by Monte Carlo simulations 

It has also been pointed out that there exists a deep 
connection between the multifractal spectrum for the 
critical wave function of a Dirac fermion in a random 
magnetic field and the spectrum of primary fields in non- 
unitary conformal field theories with vanishing Virasoro 
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conformal charge jl3|,[l4| on the one hand, and Liouville 
field theory [[L5| on the other. It is then tempting to 
speculate that there must be a counterpart to the freez- 
ing transition characterizing GREM in these quantum 
field theories. 

In the present paper we show that the exact results of 
Ref. Q on the multifractal spectrum and, in particular, 
the existence of a phase transition, can be obtained com- 
pletely within the framework of the Gaussian field the- 
ory in two-dimensional Euclidean space. The key idea 
of our approach is the introduction of a quantity fl(E) 
that counts the number of points where the critical wave 
function has amplitude (or height) \& oc e~ E / 2 in a disk 
of radius L (suitably discretized as a lattice of spacing 
a). We show that this quantity is directly related to the 
microcanonical density of states of the GREM and which 
we will also call the density of states. The advantage of 
the approach that we present in this work is that it is 
direct and it does not rely on the use of either replicas or 
supersymmetry. Our results are essentially rigorous, ex- 
cept for a conjecture, which we believe to be true, about 
the self-averaging property of the probability distribu- 
tion of the entropy hifl(E) in the thermodynamic limit 
N — > oo. We also explain in detail why the Gaussian 
field theories in both Euclidean and ultrametric spaces 
give exactly the same results. 

This paper is organized as follows. In Sec. |lj we intro- 
duce the model and the quantities of main interest: the 
multifractal exponents r(q) of the critical wave function, 
its Legendre transform f(a) and the density of states 
Q(E). Here we draw an analogy between the computa- 
tion of the inverse participation ratio of order q of the 
wave function and the partition function for all spatial 
configurations of a static particle in a random potential 
</>(x) = — In | \& (x) | . Q(E) is the microcanonical density 
of states of this equivalent problem. In Sec. Ill we derive 
the thermodynamic properties of the equivalent system. 
Here we calculate the average density of states and use 
it to show that the equivalent problem has a phase tran- 
sition at a critical "energy" determined by the width g 
of the probability distribution of the random vector po- 
tentials. We show that the thermodynamic functions of 
the equivalent problem are exactly those of the GREM. 
In particular, we show that, in terms of the critical wave 
function, the phase transition of the GREM represents 
the onset of the regime where the probability distribu- 
tion of the wave function is undersampled in a given dis- 
cretization of the plane. In Sec. [V we use these results 



to derive the exact form of the functions f(a) and r(g). 
Our results agree completely with the analysis given in 
Ref. B. Section [y| is devoted to the conclusions. Techni- 
cal details of our calculations are given in the appendices. 



II. MODEL 

We consider the problem of a massless Dirac fermion 
moving on a plane and interacting with a static random 
magnetic field normal to the plane P,p| Jl^ , [l^ ,|6|,|l7| . In 
this model, the wave functions are localized for all en- 
ergies other than the critical energy E = 0, at which 
the wave function is multifractal |^],[l3| • This model thus 
describes a metal-insulator transition in two dimensions. 

The Dirac Hamiltonian in random and static vector 
potentials in two space dimensions is 



(2.1) 



For convenience we will set the Fermi velocity vp to unity 
from now on. This operator acts on the space of nor- 
malizable (in a finite area) two-c omp onent spinor states 
^/q,(x), with a — 1,2. In Eq. (2.1) er denotes a two- 
component vector of two 2x2 Pauli matrices, which we 
take to be a\ and oi respectively, with fj, = 1, 2 being the 
two orthogonal directions on the plane. The probability 
distribution of the random vector potential A(x) will be 
specified below. In principle, other sorts of random fields, 
such as random mass and random chemical potential, are 
also allowed in a general situation Q . 

If only random vector potentials are allowed, E = is 
an exact eigenenergy for all realizations of the disorder 
and the corresponding wave functions can be determined 
exactly [[l8). Indeed, let be a two-component spinor 
with energy E = 0, i.e., H^fa = 0. Then by means 
of a combination of chiral and gauge transformations, 
parametrized by the fields 0(x) and x( x )j respectively, 



*o(x) 



«x( x ) 



(2.2) 



the eigenvalue equation reduces to the requirement that 
the two component spinor 77 (x) satisfies the nonrandom 
Dirac equation 



«c M 5 M ?7(x) = 0. 



(2.3) 



Here, the chiral "angle" 0(x) and the g aug e transforma- 
tion x( x ) must be chosen to solve Eq. ( |2.5| ) below. 

In this work we will only be interested in the mul- 
tifractal properties of the amplitudes of the wave func- 
tions. These properties involve only the magnitude of 
the wave functions and are independent of their phases, 
and hence are gauge invariant properties. Howe ver, it 
should be stressed that the decomposition of Eq. (2.5) is 
only valid if the total magnetic flux threading the disk 
always vanishes, an assumption that we will implement 
by an appropriate choice for the probability distribution 
of the vector potential A(x) |lS(| . The remaining degrees 
of freedom carried by the spinor ry(x) then span a two- 
dimensional Hilbert space. Here we will choose the spinor 
?7(x) = (1, 0) for convenience. It is worth noting that any 
choice of spinor breaks the chiral symmetry generated by 
<73 . This procedure makes sense if we think of switching 
on an average uniform magnetic field and then taking 
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it to zero. Indeed, a magnetic field selects a state with 
unique chirality determined by the sign of the magnetic 
field. In the context of the Chalker-Coddington model 
the choice of spinor is thus the equivalent of the choice of 
the chirality of the edge current for a system on an open 
geometry such as a disk. 

With the above considerations we write the E = 
wave function as 



(2.4) 



and drop the constant spinor rj(x) = (1,0) altogether. 
The random vector potential and magnetic field are given 
by: 



A v (-x) = E„ p d p (f>{x) + c?„x( x )> 
B(x) = -V 2 0(x). 



(2.5) 
(2.6) 



Notice that the gauge degrees of freedom enter through 
X and not </>, and any phase in the wave function can 
be eliminated by a gauge transformation. Finally, we as- 
sume a Gaussian distribution of magnetic fields p(J as 
follows: 



P[0(x)] (X 



-i / d 2 x(V«;Hx)) 2 



(2.7) 



where g is the width of the probability distribution and 
plays the role of a coupling constant in this problem. One 
verifies that, in the thermodynamic limit, the uniform 
magnetic field does indeed vanish for all realizations of 
the disorder. As a technical point, it is understood that, 
among all fields 0(x) that differ from each other by a 
uniform value, only one representative is counted in the 
disorder average over <p [for example, of all possible uni- 
form fields, only the field 0(x) = is counted]. In this 
way, the ambiguity in the many to one relation between 
</>(x) and A(x) in Eq. ( |2.5| ) is removed (see Appendix |a|). 

The multifractal nature of the wave function can be 
probed through the moments of the probabilities p x ob- 
tained from the normalized wave function W(x). As 
we anticipated above, we shall use a lattice regulariza- 
tion with a microscopic cutoff distance a (box size), and 
macroscopic system size L. There are thus L 2 /a 2 sites. 
The moments of p x are the inverse participation ratios 
V(q,a/L): 



X 
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J2 e -2<?0(x) 
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The multifractal exponents r(q) for a given wave function 
arc then defined by fBEft 



r(q) = D(q)(q-l) 



lim 



hiV(q,a/L) 
Hi) ' 



(2.9) 



We will show later that r(q) is self-averaging, i.e., that 
for any realization of the disorder r{g) — (r(q)), where 
(■ ■ •) denotes the avera ge o ver the distribution of random 
magnetic fields of Eq.(2.7). 

An equivalent way of describing the multifractal prop- 
erties of a given wave function is through the scaling ex- 
ponents a x , defined by @-p5[ 



Px 



(ir 



(2.10) 



The number of lattice points at which the exponent a x 
takes values between a' and a' + da' defines the functions 
p and /: 



doi 



-/(«') 



(2.11) 



As is well known |2jj2^], the two sets of exponents r(q) 
and f(a) are related by a Legendre transformation: 



Q _ dr(q) 
dq 

f(a) = aq 



r(q). 



(2.12) 
(2.13) 



It is possible in general to map the quantities that de- 
scribe a given multifractal wave function into thermody- 
namic quantities [^6| : q maps into an inverse tempera- 
ture /3, r(g) into a free energy /o(/3), ot into an internal 
energy e, and f(a) into an entropy s(e). In our case, 
however, the equivalence is more evident due to the par- 
ticular form of the wave function of Eq. (|2.4|). We map 
our problem into the statistical mechanics of a single par- 
ticle in a lattice of spacing a and size L with a random site 
potential V(x) = 20(x), in the static limit (i.e., for the 
hopping matrix element equal to zero). For this model, 
the random canonical partition function for a particular 
realization of the disorder reads: 



(2.14) 



where the role of the random energies is played by the 
values that the disorder potential V^x) = 20(x) takes. 
For each disorder realization, the free energy of the sys- 
tem is 



HP) 



\nZ(/3). 



(2.15) 



In this system, the number of thermodynamic degrees of 
freedom is 



3 



.V = Inl- 
et 



and the number of energy levels is 

2 



(2.16) 



(2.17) 



This is similar to a system of N spins with S = h , where 
one has 2 N states. Thus, the intensive free energy is 
defined as 



MP) 



= F(P) 
N ' 



(2.18) 



[The subscript is introduced to avoid confusion with 
the spectrum f{a)] 

In this problem it turns out that it is also useful to 
count states directly, which naturally leads us to define 
a microcanonical partition function (or density of states) 
for each random field configuration: 



(1(E) = ^ <V [E - V(x)] 



(2.19) 



Here S\y(E) counts the number of states in a region of 
width W around E. One can choose for it, e.g., either 
a top hat (i.e., a product of step functions) 5w(E) = 
6{E - W/2) 8(W/2 — E), or a smooth function such as 



S W (E) 



(2.20) 



The microcanonical and canonical partition functions are 
then related by Laplace transformation: 

Z{J3) = J ^ (1(E) e^ E . (2.21) 

Here, as is usual in statistical mechanics, W has to be 
taken as small as possible but still larger than the av- 
erage level spacing A. As we show later, in the present 
case this translates into the condition: 



1» W» A = 4 A /f ^ 

2?r (L\ 



(2.22) 



From the microcanonical partition function one can ob- 
tain the total entropy: 



S(E)=ln(l(E), 



(2.23) 



and again it will be convenient to define the entropy and 
energy per thermodynamic degree of freedom: 



S/N, 



E/N. 



(2.24) 



Having established the thermodynamic dictionary, the 
normalization factor for the wave function and the inverse 
participation ratios V(q, a/L) can be written in terms of 
the partition function: 



*(x) 



(x) 



V(q, 



L' 



VW) 

Z(q) 

z(i) q - 



(2.25) 
(2.26) 



In turn, the function r(q) can be simply expressed in 
terms of the free energy at the inverse temperatures (3 — q 
and (3 = 1: 



r{q) 



HQ) ~ F(l) 

Hi) 

2q \ lim [/o(«) - /o(l)] • 



-q lim 



(2.27) 



Similarly, from the definition of a x [see Eq. (2.10)] and 
the value of the wave function [Eq. (|H)], we see that 

7(x 



/ (i; 



(2.28) 



Thus, comparing the definitions of f(a) [Eq. ( 2.11 )] and 
the entropy, one verifies that 



o = 2[e-/ (l)], 
f(a) = 2s(e). 



(2.29) 



By using the relation between r(q) and fo{P), and be- 
tween /(a) and s(e) one recovers the well known result 
that the Legendre transfo rma tion in the language of mul- 
tifractal exponents [Eq. ( [2.131 )1 is equivalent to the ther- 
modynamic relation: 



fo(P) 



£(e) 
/3 ' 



(2.30) 



III. THERMODYNAMICS OF THE MODEL 

In this section, we are going to show that the proba- 
bility to find energy levels outside a particular window 
of energy [— e c , e c ] vanishes in the thermodynamic limit. 
This is the most important result of our work. 

The key step in our argument is our estimate for the 
disorder average of the microcanonical partition function, 
which relies crucially on the choice of a distribution for 
the disorder with a variance depending logarithmically 
on spatial separation. 

The simplest way to study the thermodynamic prop- 
erties of the system is to obtain the disorder average of 
the microcanonical partition function: 

((1(E)) = [vV(x)P[V(x)}^25 w [E-V(x)} 

J X 

P0(x)P[0(x)]^<V[£-20(x)] 
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T ^ ex P{-^ 2 / 2 [^ 2 + 4G(x,x)]} 



x /W 2 +4G(x,x) 



■L\ 2 e- E2 ^^)} 
W I — 
a 



7r V a / 



(3.2) 



Here, G(x, y) is the Green function | E | for the field the- 
ory defined by the action of Eq. fl2.7| ), with a short dis- 
tance cutoff length given by a. Its form in the short 
distance limit is 



G(x,y) = -|-ln 



|x-y r 



v- 



(3.3) 



On the last line of Eq. (f|j) , we have neglected the width 
W as compared to terms of order In (— ). 

Alternatively, in terms of intensive quantities, the dis- 
order averaged number of states in an energy interval of 
width w — W/N around e = E/N is: 



(0( 



N jvfi 

W , I e V 4(3/2-), 

2 ff 



(3.4) 



This means that in the thermodynamic limit (N — > oo), 
(0(e)) goes to zero for |e| > 2-^/g~/2n and diverges ex- 
ponentially with the number of degrees of freedom for 
|e| < 2yp/27T. This behavior for the average number of 
states (0(e)) indicates that some sort of transition should 
occur at the critical energy |e| = e c where 

(3.5) 

To show that, indeed, this is a phase transition, we must 
consider not the average (0(e)), but the number of states 
0(e) for a given realization. Indeed, as shown in Ap- 
pendix g, n(E) is a random variable with very strong 
fluctuations. Hence, we cannot rely, a priori, on any 
single moment of fl(E) to identify a phase transition. 
Rather, we must use stronger probabilistic arguments to 
show that e c is indeed the critical energy of interest. We 
do so below by considering the two regimes |e| > e c and 
|e| < e c separately. 



A. |e| > e c 

In this region the thermodynamic limit of the micro- 
canonical partition function is (0(e)) — > 0. Alternatively, 
we may say that for a large enough system size N, the 
average number of states (0(e)) <C 1 in this region. If we 
naively try to obtain the entropy as ln(0(e)), we would 
find that it becomes negative for |e| > e c . This is not 
correct; the entropy must be defined for each and every 



realization of the disorder. That means we must focus 
on S(e) = lnil(e) for separate realizations. 

We are going to show now that the probability of find- 
ing any state at all with |e| > e c vanishes in the thermo- 
dynamic limit. To do that, we define the random variable 
0>(e) that counts the number of states with energies e' 
> e: 



0>(e) = 



"de>^ + 



de 



,0(e') 



w 



(3.6) 



Due to the fact that Jl> (e) is either positive or zero, we 
can bound the probability P{0>(e) > 1} of finding at 
least one state with |e'| > e by the average of 0>(e) [E7|: 



P{0>(e) > 1} < (0>(e)). 



(3.7) 



Equation ( |3.7| ) is very general since it applies to any ran- 
dom microcanonical partition function. However, this 
inequality becomes very powerful when combined with 
our estimate Eq. (3.4) for (51(e)). Indeed, with the help 
of 



(fi>(e))=-/ de'(0(e')) 



(3.8) 



^W(e/e c ) 



-N(e/e c ) 2 -l] 



\NeJ 



and our upper bound Eq. (3.7), we see that P{0>(e c ) > 
1} vanishes in the thermodynamic limit N — > oo. We 
conclude that, for a given realization, the energy levels 
will fall within the interval |e| < e c with probability 1 in 
the thermodynamic limit. 

It is now cl ear h ow to estimate the average level spac- 
ing A in Eq. ( 2.22 ). We simply multiply 2e c by the ratio 
of the number TV of thermodynamic degrees of freedom 
to the number e N of energy levels: 



A = 



2iVe, 



,N 



= 4 




(3.9) 



Finally, we can bound from above the probability 
for the density of states 0(e) to be nonzero. For a 
given realization, f2(e) counts states in an energy interval 
w = W/N around e. Therefore 0(e) must be a positive 
integer |ps| , in contrast with the average (0(e)). In the 
same way as before, the probability of having a nonzero 
microcanonical partition function is bounded from above 
by a number that goes to zero in the thermodynamic 
limit WA- 



P{0(e) > 0} = P{0(e) > 1} < (0(e)). 



(3.10) 



The entropy is therefore not defined for |e| > e c . 

One can get some intuition for these results with one 
word: under sampling. Consider f2(e) for an individual 
realization as the histogram of the number of states per 
energy interval. One can think of 0(e) as the product of 
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two factors: one is a constant that counts the total num- 
ber of states in the system, and the other is the probabil- 
ity (normalized to 1) for one state to fall within a given 
energy interval. In the present case the first factor has 
the value e N and the second is Gaussian with a width 
that grows linearly with N. In other words, the number 
of data points we have is e N , which is not enough to sam- 
ple the tails of the Gaussian distribution. If the number 
of data points were e N1 , with 7 > 1 , by taking N large 
enough we could sample the whole distribution, and all 
bins in the histogram would contain a positive number 
of points. 



B. \e\ < e c 

In this case we can define the entropy for a given re- 
alization S(e) = In 17(e) [recall Eq. fl2.23| )]. We will now 
show that the entropy, for all realizations of the disorder, 
has a common upper bound that scales linearly with sys- 
tem size. Furthermore, we will also show that the prob- 
ability to find a realization with entropy less or equal 
than ln(f2(e)) is equal to one in the N — > 00 limit. The 
arguments go as follows. 



According to our definition Eq. ( 2.19 ), 51(e) counts the 
number of energy levels in some energy window. That 
number cannot be larger than the total number of energy 
levels e N = L 2 /a 2 . Consequently, lnf2(e) is bounded 
from above by N. Incidentally, no such bound holds for 
ha.Z(j3) for individual realizations of the disorder. 

Next, we introduce the probability to find the intensive 
entropy s = (1/N) lnfi(e) in a given interval (si, S2) by 



P/v(si < s < s 2 ) = / dfj, N (s) 



(3.11) 



Here, d^Ar(s) is the measure of the entropy for N degrees 
of freedom. We will make the following assumptions: 

1. The probability measure d/ijv(s) has a well defined 
thermodynamic limit: 

lim / dfi N (s) = / dfioc(s) (3.12) 

for all s% < S2- 

2. There exist two positive constants A and B such 
that 



+00 



(n(e)) = / d/J, N (s)e = vNe 



AN+B 



(3.13) 



for N sufficiently large. 



Notice that we have proved the validity of the second 
hypo the sis and that A and B can both be read off from 
Eq. (3.4). In particular, 



(3.14) 



We now choose s+ > A. By assumption, s_|_ > and 



(fi(e))> / dn N (s)e 



r+oc 

> / d^ N (s)e 1 

Js + 

= e Ns +P N (s+ < s). 



(3.15) 



We have thus established the upper bound, 

Pn(s+ < s) < VNe^ s ^ N+B , (3.16) 

for all s + > A and for all sufficiently large N. Since the 
thermodynamic limit is assumed to be well defined, we 
conclude that 



Poo(A < s) = 0. 



(3.17) 



In other words, the probability to find the entropy density 
(l/iV)lnfi(e) > (1/N) ln(0(e)) vanishes in the thermo- 
dynamic limit N — > 00. 

If we could prove that the probability to find an in- 
tensive entropy lower than (1/N) ln(f2(e)) also vanishes 
in the thermodynamic limit, we would have shown that 
the entropy density is indeed self-averaging (i.e., al- 
most surely takes one and only one value), and given 
by limjv->oo(l/-^0 ln(f2(e)). However, without further 
knowledge of the probability distribution for the entropy, 
it is not possible to prove the existence of a similar lower 
bound. In fact it is possible to c onstru ct a sequence of 
measures which obey Eqs. (3.12) and ( 3.13j ) but is not 
self- averaging in the thermodynamic limit. We neverthe- 
less expect the entropy density to be self-averaging for 
two reasons. First and in addition to our upper bound, 
it is possible to estimate the ratio of the n-th moment 
of the partition function Z(f3) to its mean raised to the 
power n. This is done in Appendix |b[ This estimate 
suggests that the field theory computations of the multi- 
fractal dimensions (with replicas, supersymmetry and Li- 
ouville field theory) are reliable g] in the regime |e| < e c . 
Second, the entropy density for any GREM is known to 
be self-averaging in the thermodynamic limit [p|-|l2| and 
in view of the close connection between GREM and our 
problem we expect this to be also true here. 

Thus, although we do not have a rigorous proof, 
we conjecture that (1/7V) lnf2(e) is a random variable 
with vanishing variance in the thermodynamic limit 
N — * oo, i.e., self-averaging, and which converges to 
limjv— >.oo(l/A0 ln(f2(e)) . This conjecture is supported 
by the fact that, in addition to our bound, it is in this 
regime that the field theory computations of the mul- 
tifractal dimensions (with replicas, supersymmetry and 
Liouville field theory) are reliable. In other words, the 
entropy should be self- averaging in the thermodynamic 
limit and it is given by 



G 



S(e) = In fl(e) = (In 0(e)) = ln(ft(e)) 



(3.18) 
(3.19) 



for any realization of disorder. Hence, for this range of 
energies |e| < e c we may commute the order of In and 
(•••), i.e., the "quenched" and "annealed" entropies co- 
incide [El. 



C. Temperature and free energy 

Up to this point we have proven that the thermody- 
namics of this model is identical to that of the random 
energy model || for |e| > e c and argued that this is 
also so for |e| < e c . In order to obtain the free energy 
as a function of temperature, we follow Derrida ||. In 
the region that contains states (|e| < e c ), we obtain the 
temperature from 



ae 



2e' 



(3.20) 



In this region the temperature will be in the range 
|T| > e c /2, and the free energy, computed from the rela- 
tion dfo/dT = — s, is given by 



/ (T) = -T-^. 



(3.21) 



In the limit when T — > T c = e c /2, we obtain e — ■+ — e c , 
s — ► and /q — * — e c . Below this temperature, the sys- 
tem cannot lower its energy because there are no states 
for e < — e c . The system remains frozen at e = — e c , 
with entropy s — and free energy fo(T) = — e c . In 
other words, for |T| lower than the freezing temperature 
T c , the energy, the entropy and the free energy remain at 
their accumulation points. In summary, 



- 1 



ol/3| 
9c ' 



if l/3| < 9c 
if |0| > q c 



(3.22) 



where q c = y^2ir/g. 

The behavior in the low temperature regime becomes 
clear if we consider the \J3\ — ► oo limit. For any configu- 
ration and for a given value of L/a, there is an absolute 
minimum -Emin and an absolute maximum £ max for the 
values of the V(yi). For large enough \(3\ these extreme 
values will dominate the partition function, and one ob- 
tains, in the case of positive sign for (3: 



-0E„ 



This means that 



lim 



fo(P) 



lim 

N— >oo 



E n 



N 



(3.23) 



(3.24) 



because we know that for a large enough system 
Emi n /N — — e c with probability one. This implies that 
fo{P) is self-averaging and equal to — e c = —2/q c , as 
shown above. 

Notice, however, that in the low temperature regime 
quenched and annealed values for the free energy are dif- 
ferent: (In Z(0) ) ln(Z(/3)). Indeed, combining Eqs. 
(f[|l]) and (|J) one obtains for all (3 the so-called log- 



normal spectrum 



\n{Z{(3)) =N 1 + 



i2 



(3.25) 



which shows no obvious sign of the phase transition. 

On the other hand, in the high-temperature regime, 
the quenched and annealed free energies (lnZ(/3)) and 
In (Z(f3)) do coincide. This can be understood by com- 
puting lnZ(/5) for one particular realization of the disor- 
der [see Eq. ( p72lt )]: 



Z{(3) = N 



de e -^-s(e)]. 



(3.26) 



If one accepts our conjecture that the entropy density is 
self- averaging in the thermodynamic limit, then the sad- 
dle point approximation on the integrand in Eq. (3.26) 
probes the parabolic part of the entropy, and one ob- 
tains for InZ(P) the result of Eq. (3.25). In other words, 
not only is \nZ([3) self-averaging but in this case In and 
(• • •) commute. In the language of multifractality, this is 
the region in which the parabolic approximation is exact. 
Besides, since at \J3\ = q c the minimum of /3e — s(e) falls 
at the edge of the populated region, there is a singularity 
in the derivative of InZ(P), and one recovers the phase 
transition point. 



D. Equivalence with random energy models 

We have calculated the thermodynamic functions of 
our model and they are identical to those of a special 
generalized random energy model || . 

It is interesting to see how this comes about. All the 
results we obtai ned really follow from our calculation of 
(Q(E)) in Eq. ( |3.l| ). The analogous calculation for the 
generalized random energy model of Derrida and Spohn 
]io| , p| gives the same result, except that positions x in 
the lattice are replaced by directed paths V in a Cayley 
tree, and the value G(x, x) = (g/2n) ln( — ) is replaced 
by G{V,V) = (g t /\nK)\nd{V,V). Then,°the two aver- 
age partition functions must correspond if the parame- 
ters are chosen properly. Note that the requirement that 
G(V, V) depends logarithmically on the ultrametric dis- 
tance dtV, V) on the tree uniquely defines the GREM in 
the thermodynamic limit. 

One can take this analysis further by studying mo- 
ments of the microcanonical partition function. Al- 
though at first sight the expressions for the two models 
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appear different, if one changes variables into logarithms 
of the distances, one can actually see that the expressions 
for the GREM actually give the Riemann sums that cor- 
respond to the integrals in the case of our model. 



IV. CONSEQUENCES FOR THE 
MULTIFRACTAL SPECTRUM 

Now that we have the values of the thermodynamic 
functions s(e) and /o(/3), we can transla te th em int o the 
language of multifractality, using Eqs. (HI) and &2$) 
as our dictionary. 

The results are as follows: starting by the spectral 
weight function /(a), we find that it is defined only in 



the interval d_ < a < d + (corresponding to the entropy 



being defined only in the interval — e c 
has the value 



/(«) = 8 



(d+ — a) (a — gL) 
{d+-d_Y ■ 



< e < e c ), and 



(4.1) 



However, the values of g?_ and g?+ will change with the 
strength of the disorder. 

Depending on the strength of the disorder, there are 
two regimes: in the weak disorder regime, which cor- 
responds to g < 2it, the quenched and annealed aver- 
ages for the logarithm of the wave function normalization 
factor Z{\) are coincident, while in the strong disorder 
regime, which corresponds to g > 2ir, they are not equal 
anymore. 

In the weak disorder regime, the extremal dimensions 
d- and d+ are both positive: 



d± = 2 



and r(q) has the form 
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r(q) 




if M < <7o 
if M > <7c- 



(4.2) 



(4.3) 



On the other hand, in the strong disorder regime, the 
lower extremal dimension is zero: 



(4.4) 




-(g) = 



tM-\q\) 



if \q\ < q c 
if \q\ > q c 



(4.5) 



Notice that in this regime we find r(q) = for q > q c , 
meaning that for all integer moments the inverse partic- 
ipation ratio does not scale with system size. This is 
usually interpreted as characteristic of a localized wave 
function. 



V. CONCLUSIONS 



We have calculated the multifractal scaling exponents 
of the critical wave function for two dimensional Dirac 
fcrmions in the presence of a random magnetic field. 
There is a transition in the multifractal spectrum, which 
is interpreted as a freezing transition common to glassy 
systems. This freezing transition is a rigorous result of 
the present work and had been previously conjectured on 
the basis of a comparison with GREM || . 

We have proven that a previously proposed mapping 
M between these multifractal properties and the thermo- 
dynamics of a generalized random energy model describ- 
ing directed polymers on a Cayley tree is indeed 
exact in the glassy regime. Our proof generalizes entropy 
considerations on GREM (Ref. 0) to a two dimensional 
Gaussian field theory. 

Derrida has also shown that a direct computation of 
the quenched free energy was possible on GREM. This 
suggests that the same could be done on the field the- 
ory. In fact, it can be shown that the generalization of 
Derrida's calculation for GREM naturally leads to esti- 
mating the partition function in Liouville field theory. 
We thus believe that there exists a counterpart to the 
freezing transition of GREM in Liouville field theory. It 
is an interesting question to probe this issue further. 

Another open issue is the fate of replica symmetry if 
the replica approach is used to calculate the multifractal 
scaling exponents. Indeed, it is known that the freezing 
transition in GREM is associated to replica symmetry 
breaking J||5(|. It would be interesting to see how this 
replica symmetry breaking manifests itself in a replicated 
version of our Gaussian field theory. 



VI. ACKNOWLEDGMENTS 



We are indebted to J. Chalker, L. Balents, 
A. W. W. Ludwig, M. P. A. Fisher, and X.-G. Wen for 
many useful comments and penetrating questions. This 
work was completed while three of us (CCC, CM and 
EF) were participants at the Program on Quantum Field 
Theory in Low Dimensions at the Institute for Theoreti- 
cal Physics of the University of California at Santa Bar- 
bara. We are very grateful to the organizers of the pro- 
gram and particularly to Professor Jim Hartle, Director 
of ITP, for his warm hospitality. This work was sup- 
ported in part by the NSF through the Grants Nos. NSF 
DMR94-24511 and NSF DMR-89-20538 at the Univer- 
sity of Illinois at Urbana-Champaign (H.C., C.C.C., E.F. 
and P.G.), NSF DMR94-11574 at Massachusetts Insti- 
tute of Technology (CM.), and by NSF PHY94-07194 at 
the Institute for Theoretical Physics of the University of 
California at Santa Barbara. 



8 



APPENDIX A: SYMMETRY AND GREEN 
FUNCTIONS 



The probability distribution for the disorder, Eq. (2.7), 
allows for an exact symmetry under a constant real shift 
of the field configuration </>(x) : 



0(x) 
1&(x) 



/>« + C, 
r^(x), 



which leaves both P[0(x)] and \l/(x) unchanged, 
ever, neither Z(q) nor Q(E) are invariant: 

n(E) -> o(E-2C), 

Z(«) -> e- 2 ^Z(q), 
(lnZ(g))-»-29C + (InZ(g)). 



(Al) 
(A2) 

How- 



(A3) 
(A4) 
(A5) 



This means that, although (In Z(q)} — q(\nZ(l)) is well 
defined and, in principle, one can compute r(q), in a 
naive calculation (\nZ(q)) would be ill defined. 

To perform an actual calculation it is convenient to 
break this symmetry in a controlled way. The simplest 
approach is to add a mass term m = l/L to the action. 
This penalizes configurations for which 0^0. Another 
possibility is to impose Dirichlet boundary conditions on 
4>(x) ■ Although Green functions for </>(x) in each of these 
cases will be different, the only value we need for our 
purposes is their short-distance limit, which is the same 
in all cases, namely, 



G(x,y) 



9_ 

47T 



In 



L 2 



(A6) 



Here, Gki is a shorthand notation for the Green function 
in Eq. ( |A6| ) with arguments x^ and xj. We notice that 
the integrand on the right-hand side of Eq. ( p32| ) does not 
depend on Xi, . . . ,x„ for n = 1 but does for n > 1. For 
higher moments than n = 1, the statistical correlations 
encoded by Gki ^ for k ^ I, imply that, in a finite sys- 
tem, (Q n (E)) ^ (Cl(E)) n . In the thermodynamic limit, 
the difference between (fl n (E)) and (ft(E)) n disappears 
if the statistical correlations are short range, as is the 
case in the original random energy model (REM) of Der- 
rida, where all energies are assumed to be identically and 
independently distributed random variables j^]. This is 
not so, however, if Gki encodes long range statistical cor- 
relations as is the case here. For any moment n > 1, 
there will be a critical energy density e c (n) < e c (n — 1) 
such that the ratio (Q n (E)) / (Q(E)) n diverges in the ther- 
modynamic limit for |e| > e c (n). In other words, f2(e) 
is broadly distributed in the thermodynamic limit. Our 
claim is that it is the limit e c = lining e c (n) that con- 
trols the freezing transition, and not the naive replica 
limit lim n _>o s c {n) = oo. 

Instead of calculating the sequence e c (n) explicitly, we 
calculate the ratio 



Rn(P) 



(Z n ((3)) 

W)>"' 



(B3) 



We find that 



APPENDIX B: MOMENTS OF n(E) AND Z(f3) 

In this appendix, we estimate moments of the density 
of states fl(E) and of the part ition function Z(f3) which, 
we recall, are related by Eq. (2.21). Such moments are 
needed to decide if quenched and annealed averages are 
equal. We begin with fl(E). Let n be an integer larger 
than one. By definition, 



Urn / D[^)]P[^)]J2- 

xi 



dXi 



(Bl) 

d\ n 



x exp 



W 



k=i 



k=i 



k=i 



We assume that all sums and integrals can be freely inter- 
changed. Averaging over disorder is a Gaussian integral 
yielding 

™>=feE-E/t-/^ (B2) 



xexp 



2 Z-^ 

fej=i 



\k(W 2 S kl + 4Gm)\i + iEJ2 



fc=i 



L 2 "' L 2 

c n (gP 2 ) 



Xj x 7 - 



L 



a X 2(«-i)(i-^#) 
L 



(B4) 



In the thermodynamic limit, the right hand side is a 
finite number for n < (2ir)/(g[3 2 ) (assuming n > 1). 
In this case Z n (f3) fluctuates weakly. But for n > 
(2tt)/( 5 /3 2 ), R n (P) diverges, and thus Z n {f3) fluctuates 
strongly. The re are important consequences that follow 
from Eq. (| 



1. There exists a sequence of critical (3 n given by 



Pn — ' 

ng n 



(B5) 



below which ln(Z n (f3))/N = ln{Z(/3)) n /N in the 
thermodynamic limit N — » oo. Remarkably, the 
same sequence of critical moments is shared by the 
random energy models studied in Refs. p|-|Tl|j, 
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2. Caution is needed when using the replica trick 

(Z"(/3)> - 1 



(\nZ{f3)) 



lim ■ 



(B6) 



Indeed, (Z n ((3)) is not an analytic function of n 
due to singularities at (3 = 2-k/ (3 2 g, and n = 1, and 
caution must be used when using the replica trick 
to calculate \nZ((3). 

For GREM, (3 C indicates a phase transition between 
the regime (3 < f3 c and the regime (3 > (3 C . For all 
GREM, the quenched and annealed free energy are equal 
if (3 < /3 C . However, quenched and annealed free energy 
are not equal if (3 > (3 C , and need not obey the same 
functional dependency on f3 for different GREM in this 
regime of temperatures. There are essentially two Gaus- 
sian GREM's whose quenched free energy obey the same 
functional dependency on (3: the REM with uncorrelated 
and the GREM with logarithmic correlated 
TT| . Both undergo a sharp freezing transi- 
tion characterized by a discontinuous one step specific 
heat (as opposed to continuous). Since we have proven 
that the same freezing transition characterizes our mul- 
tifractal scaling exponents, we conclude that they must 
be self-averaging and given by an annealed average below 
the critical moment q c . 
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